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Abstract 

Band theory predicts both K3C60 and K4C60 to be metals; various experi- 
mental probes show that while K3C60 is indeed metallic, K4C60 appears to be 
insulating. The standard view of this apparent failure of the single-particle 
picture is that electron correlation is predominant. We describe an alter- 
native scenario, motivated on theoretical grounds, which invokes a spin- or 
charge-density-wave state to explain the observed insulating behavior. 

(To appear in Physica B) 
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I. INTRODUCTION 



Since the discovery of superconductivity in the potassium-intercalated fulleride K 3 C 60 , 
considerable theoretical attention has focused on the electronic properties of the A u Cqo 
family (A=alkali), much of it assuming the validity of a single-particle band description. A 
growing body of evidence suggests that this assumption may rest on rather shaky ground, at 
least for some members of the family. A striking example is the contrast between K3C60 and 
K^Cqq. According to first-principles band calculations, both are metallic, with comparable 
Fermi-level densities of states. Photoemission gives a very different picture, however: for 
K3C60, strong Fermi-level emission is observed, while for K4C60 the spectrum is insulating 
Q. Additional evidence is found in NMR, which also shows that K3C60 is metallic and K4C60 
is insulating || . Several workers have argued that strong electron-electron correlation splits 
the conduction band (CB) into a filled lower Hubbard band and an empty upper Hubbard 
band 

In this paper, we propose a different scenario to explain the insulating nature of the A^Cqo 
fullerides. We begin by discussing the results of first-principles band-structure calculations 
for K4C60, then explore possible broken-symmetry ground states such as a charge-density 
wave (CDW). We emphasize that our results do not in any sense prove that a CDW state is 
the thermodynamic ground state; rather, we make a plausibility argument and then describe 
model-Hamiltonian calculations which may give further insight. 

II. BAND STRUCTURE OF K 4 C 60 

We use the local-density approximation (LDA) to density-functional theory. Bloch basis 
functions are linear combinations of occupied and unoccupied atomic orbitals for potassium 
and carbon; these are in turn expanded on a set of gaussian functions. This provides a 
compact basis, and allows us to perform accurate all-electron calculations — without the 
need for pseudopotentials — in which core, valence, and conduction states are treated on 
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equal footing. The potential and charge density are completely general and without shape 
approximation. We have previously used identical methods to study K 6 C 60 |§, K 3 C 60 0, and 
B&sCqo @; further details concerning numerical methods may be found in these references, 
and preliminary work on ^Ceo may be found in Ref. |[7||. 

The K4C60 geometry used in our calculation is taken from the x-ray refinement of Fleming 
et al. J8|, who found a body-centered tetragonal (bet) lattice for the A^C^o family, with lattice 
constant a=11.886 A and c/a=0.906 for y4=potassium. The potassium positions were refined 
to (0.22,0. 50, 0)a and symmetry-related sites. Fleming et al. observed that because Ceo lacks 
a four-fold symmetry axis, a tetragonal lattice (a = b) implies that the Ceo molecules are 
either orientationally disordered or are spinning rapidly; Zhou and Cox have shown that 
a merohedrally disordered model fits the data well ||. Rabe has proposed an alternative 



structural model fllP] , consisting of an orientationally ordered crystal with a screw-axis 
symmetry that leads naturally to a = b. Here, we make the standard simplification of 
taking the K 4 C 60 structure to be orientationally ordered, with one C 60 molecule per cell. 
Work currently in progress, which goes beyond this simplification, will be described in 
Sec. 0. 

In Fig. [I] we show the band structure for orientationally ordered K 4 C6o- Since the bet 
conventional cell is a slightly distorted version of the body-centered-cubic (bec) cell of K 6 Ceo, 
we show the bands along the bet route that corresponds most closely to the bec route appear- 
ing in Fig. 2 of Ref. [§J. Indeed, the K 4 C 60 bands bear a close resemblance to those of K 6 C 60 , 
with an important difference. The point group of orientationally ordered A^Cqq (generated 
from three Cartesian mirror planes) has only one- dimensional irreducible representations, so 
that any band degeneracies must be accidental. This has the effect of smearing somewhat 
the spectral density throughout the CB, resulting in a DOS with few sharp features. The 
Fermi-level DOS N(ep)=18.3 states/eV-cell (both spins) is ~40% larger than for K3C60, 
strongly at variance with the NMR and photoemission results discussed above. 

The Fermi surface (FS) has 3 sheets, shown in Fig. ^. Two form closed hole pockets 
centered at T. The third forms quasi-two-dimensional sheets, large parts of which are nearly 



flat over much of the zone and roughly normal to the c-axis. To the extent that these 
sheets are truly flat, they represent an unusual form of FS nesting, with a nesting vector 
Q = Qz. It is well known that FS nesting can lead to electronic instabilities, by enhancing 
the susceptibility x(Q) f° r the formation of a charge- or spin- density- wave (SDW) ground 
state. The simplest illustration of this occurs in a linear chain of atoms equally spaced by 
a, for which a half-filled cosine band has a nesting vector equal to half the zone, Q = it /a. 
The Peierls mechanism provides a route to lowering the total energy, by chain dimerization 
(that is, by forming a commensurate CDW with period 2a) and the opening of a gap. 

To quantify the degree of nesting and identify posible nesting vectors in K4C60, we 
calculate 

£nm(Q) = iV _1 ^5(e nik - e F )S(e mM+Q - e F ), (1) 

k 

which measures the phase space available for any process that scatters electrons from point k 
on the nth FS sheet to k+ Q on the mth sheet. This factor also enters into the second-order 
susceptibility, d 2 E t /dhQdh_Q, of the total energy E t with respect to density variations /iq, 
making explicit the connection between nesting vectors and CDW wave vectors. In Fig. |3| 
we plot the function £33 (Q) in the Q x , Q z plane, revealing a single moderately strong nesting 
vector centered near Q = ±(0,0,7r/c). This is precisely the nesting vector required for the 
smallest possible commensurate density wave, namely, period doubling along the c-axis. In 
the following section we describe preliminary work, using a model Hamiltonian, that may 
lend further insight into this scenario. 



III. MEAN-FIELD HUBBARD MODEL 

Shortly after the discovery of superconducting K3C60, Zhang, Ogata, and Rice proposed a 
negative-f/ model for fulleride superconductivity mediated by alkali-ion optical phonons ||TT] . 



They argued that if the interaction is attractive, then superconductivity would result from 
the competition between a CDW ground state and a superconducting ground state, for which 
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the fee lattice of K3C60 frustrates the former in favor of the latter. Since this proposal, isotope 
experiments have demonstrated that alkali optical phonons are not involved in pairing. 
Moreover, considerable experimental and theoretical evidence has shown that the electron- 
electron interaction is intrinsically repulsive. Nevertheless, we discuss the negative-?/ model 
here for two reasons: (1) we will later use the same formalism to study the consequences of 
a repulsive interaction, and (2) it provides a pedagogically instructive context for a theorem 
which states that for the U < Hubbard model with a single band and nearest-neighbor 
hopping, the CDW and superconducting ground states are exactly degenerate when the 
band is half filled and the lattice is bipartite [0. For K3C60, the band is half filled but the 
lattice is not bipartite; the authors of Ref. [1]J showed numerically that the degeneracy is 



broken in favor of a superconducting ground state. For K 4 C 60 , the situation is reversed: the 
lattice is bipartite but the filling is two-thirds, again breaking the degeneracy. It is tempting 
to hypothesize that since (1) the bet lattice is bipartite, and hence no longer obviously 
frustrates the CDW state, and (2) the ratio c/a < 1 provides a preferred direction for the 
CDW wave vector, then perhaps the degeneracy might be broken in favor of the CDW state. 



To explore this possibility, we proceed along similar lines to Ref. |TT] , solving the negative 



U Hubbard model within the mean-field approximation. We begin with the Hamiltonian 

# = -* E ( c fe + h.c.) +Uj2n i ^n ii . (2) 

(*,i> * 

We assume a single band with nearest-neighbor hopping on the bet lattice; the dispersion is 
then 

/9k = —8 cos(k x a/2) cos(k y a/2) cos(k z c/2). (3) 
For the CDW state, the energy per electron is 

e cdw = ^ e <r+^[i-(vo 2 ], (4) 

^> 7V k,s=±l 4 

where 5 is the order parameter given by (rii) = 1 ±5, and ( = 4/3 is (twice) the band filling. 
The sum is over filled states and the band energy is given by 



CDW 



(tfa) 2 + (U5/2) 2 1/2 . (5) 



For the superconducting state (SS), the energy per electron is 



E ss = -±Z4 a + e F + ^[l-{2A/C)>\, (6) 
where A = (etc-,) is the superconducting order parameter, and the band energy is now 
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(Ac - e F f + (UAf 1/2 . (7) 



In both Eqs. (f|) and M), the total energy is evaluated by minimization with respect to the 
order parameter. 

We have calculated the CDW and SS energies for < \U\/t < 10; the model predicts 
the SS state to be favored throughout this range. This is strongly at odds with experiment 
(which finds no evidence for superconductivity) and with our hypothesis of a favored CDW 
state. We discuss this below. 

Although instructive, the Hamiltonian of Eq. (Q) is too crude in several important ways: 
(1) A single tight-binding band does not reproduce the nesting features of Fig. ||. (2) 
Recent work on tight-binding parameterizations of the LDA energies shows that nearest- 
neighbor interactions alone are inadequate for reproducing the LDA spectrum, and that 
two coordination shells are necessary fll3"| . (3) A variety of experimental and theoretical 
evidence has led to a consensus that U is positive, with a value of order 1 eV. To address 
these deficiencies, we have begun to study the Hamiltonian 

H = 12 { c L T ^j" c jv + h - c ) + U Y, n i\ n ni ( 8 ) 

i/j,,ju i 

where cj creates an electron on site i in orbital \i = x,y, z. The parameter t has been 
generalized here to a set of 3x3 matrices giving the amplitude for hopping between different 
orbitals on different sites (possibly between molecules in different orientations). We retain 
hopping terms between first and second neighbors. The matrix elements T ifl j u are obtained 
from the LDA eigenvalue spectrum by Fourier inversion, giving the optimal tight-binding 
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representation of the LDA dispersion. Finally, we study both positive- and negative-?/ 
solutions, i.e., both SDW and CDW states. Preliminary indications are of an enhanced sus- 
ceptibility for forming a density-wave state with Q oriented along the c-axis, corresponding 
to a value of \U\ ~ 1 eV 0. 



IV. PERIOD DOUBLING AND BAND STRUCTURE 

A complementary line of inquiry can be followed within the single-particle picture. LDA 
has previously been used to study commensurate CDW's by first assuming, say, a doubled 
unit cell and then minimizing the total energy with respect to a dimerization coordinate. 
Whether or not this procedure works in principle is the subject of debate: Overhauser has 
argued that correlation effects are important for density-wave ground states and are too 



crudely represented by LDA [T^], while Ashkenazi et al. have shown numerically that LDA 
fails to predict any dimerization even for the textbook Peierls system, polyacetylene ||15| . 

For these reasons, we focus on band-structure effects. We assume that a density wave 
in K4C60 leads to period-doubling along the c-axis (as suggested by the FS nesting), with 
an amplitude of 5%. The orientations of the Cgo molecules are assumed not to change. The 
resulting LDA band structure for this doubled cell is suggestive, but not definitive: no gap 
appears, but the extended flat FS sheet disappears, leaving only small electron and hole 
pockets. We estimate, from simple FS-area scaling arguments, that the dc conductivity 
is smaller by a factor of 5-10. While not insulating per se, the reduced conductivity may 
be manifested in NMR as a small Korringa component plus an activated component. Ad- 
ditional effects from disorder (orientational or defects) are expected to reduce further the 
measured conductivity. 



We acknowledge helpful discussions with E.J. Mele and W.E. Pickett. 
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FIGURES 

FIG. 1. Theoretical LDA band structure for (orientationally ordered) IQCeo- The Fermi level 
is the energy zero. 

FIG. 2. Fermi surface corresponding to the LDA bands of Fig. 1, in a repeated-zone represen- 
tation. 

FIG. 3. Surface plot of ^33(Q) i n the Q X ,Q Z plane. The plotting rectangle is centered on 
Q = 0, and has sides of length 3(ir/a) and 3(tt/c) in the Q x and Q z directions, respectively. The 
singularity at Q = is physically meaningless. 
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